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PREFACE 



"Underway on nuclear power." With this simple but historically 
significant message, the U.S.S. Nautilus heralded a new age in power 
generation. Powered by the first mobile nuclear reactor, she set the 
stage for dramatic developments in military and industrial fields. A 
powerful vessel, but unwieldy and undesirably large for a submarine, 
her size was largely determined by the size of the power plant com- 
ponents. One of the more significant of these components is the pres- 
surizer for the reactor primary coolant system. This tank is extremely 
large due to the conservative design necessitated Dy ignorance of the 
thermodynamic transient behavior in its pressure range. 

The objective of this thesis is to produce design data for steam 
pressurizing systems, oy electronic analog simulation of the thermo- 
dynamic transients which occur in the pressurizer vessel. 

The writers wish to thank Professor Eugene E. Drucker, and Profes- 
sor Hugo M. Martinez of the U.S. Naval Postgraduate School for their 
assistance and enthusiastic cooperation in the prosecution of this 
project. • 
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CHAPTER I 



INTRODUCTION 



The prototype (Hark 1) of the Submarine Thermal Reactor (STR) was 
the first full scale power reactor in the world to be completed and 
successfully operated. The first central station nuclear power plant 
in the United States will be the Pressurized Water Reactor (p,7K) , 
nearing completion at Shippingport , Pennsylvania. The many highly 
desirable aspects of the pressurized light water reactor cause it to be 
one of the most promising types, to date. 

A basic requirement for these light water reactors is that a very 
high pressure is maintained on the primary coolant. Pressures in the 
neighborhood of 2000 psia and higher are presently in use. High pres- 
sures permit the use of sub-cooled water at high temperatures in the 
reactor without danger of boiling. 

'When any control program other than constant average primary loop 
temperature is used, a change in the volume of the primary coolant is 
to be expected for a corresponding change in power. A volumetric tran- 
sient is consequently induced and a surge tank is required in the 
system. Since the primary coolant must be maintained at a high pressure 
the surge vessel is also used to perform this function. 

The pressurizer must be strong enough to maintain steam and water 
in equilibrium at high pressures, it must be large enough to absorb 
normal and accidental surges without permitting excessive pressures 
in the primary loop. Although thermally insulated, it must have internal 
heaters, capable of maintaining the water and steam at the saturation 
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temperature. These heaters must have additional heating capacity cap- 
able of generating steam at sufficient rate so as to prevent excessive 
pressure drop in the primary coolant during negative (out) surges. 

Specifically, in this thesis we will assume a primary loop working 
pressure of 2000 psia, normal volume surges of 4 cubic feet and accidental 
surges of 7 cubic feet as representative values. At steady state water 
and steam are maintained at saturation temperature (636 °P) in the tank, 
but only the pressure (2000 psia) is transmitted to the cooler primary 
loop via a relatively long stand-pipe. Assuming that the installed 
heaters adequately compensate for negative surges, this type surge is 
not considered a design limitation. Therefore only those surges caused 
by an increase in primary loop volume, or positive surges, will be 
considered here. 

When designing or sizing a pressurizer, limits must be imposed on 
the various properties of the fluids in the tank. The maximum allowable 
pressure change must be determined from consideration of the strength 
of the members involved. The expected change in average density and 
therefore volume of the primary coolant must be computed. The duration 
of the surge is determined from changes in power, either accidental or 
deliberate. The one remaining variable, the size of the tank must now 
be determined on the basis of the aforementioned parameters. 

The pressure change caused by a positive surge may be minimized by 
introducing a portion of the surge water into the top of the tank as a 
finely divided spray. The most conservative basis for pressurizer 
design is to assume that no „spray is employed, the steam is dry saturated, 
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and undergoes isentropic compression. Should the steam not be ary, the 
resulting final pressure would be lower than in the case of dry steam. 
The optimum and therefore minimum size required can only be positively 
determined after a complete understanding of the thermodynamic and heat 
transfer mechanisms involved. 

In this project, it is felt that valuable information is available, 
both in the fields of design limitation and in determining the charact- 
eristics of high pressure saturated steam. A Boeing Electronic Analog 
Computer with associated 4-channel Sanborn Recorder (see Fig. l) is 
utilized to simulate the pressurizer, with the above objectives in 
view. 

An important assumption which is made at the outset is that the 
driving function, specific volume (which is directly related to tank 
level) follows a cosine curve. Analysis of actual tank transients 
indicates that this approximation is a good one for a large portion of 
them. This type of function can be easily produced on the analog com- 
puter, while production of more exact functions is difficult. 

In the analysis of results, every attempt is made to utilize the 
principle of geometric similarity, to more nearly generalize these 
results . 
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BOEING ^FCLluNlC ANALOG Co^UTiid *LD SAuBum 

Figure 1 
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CHAPTER II 



DEVEL0P1AEI.T OF THEORETICAL EQUATIONS 

The initial step in the analog simulation of this problem consists 
of a description of the system in terms of theoretical equations. 

.Vhenever possible known thermodynamic relationships are used, and when 
this is not possible empirical relationships are developed, to reference 
can be found in the literature treating the following development of 
descriptive equations, and it is believed that this approach may be a 
new one. 

As an initial simplifying assumption, an equivalent system is devised 
which lends itself more readily to analysis than does the actual tank. 



PERFECT 

INSULATION 



Equivalent System 
Figure 2 
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The actual tank is well insulated, therefore the concept of perfect 
insulation for the equivalent arrangement does not introduce significant 
error. The thermodynamic system is defined as the mass of steam in the 



tank. In the actual tank tnere is installed internally a spray and 
degasifier assembly. In this analysis of positive surges, without spray, 
its only effect is to function as a heat sink. The heat sink in the 
equivalent arrangement represents the tank walls, spray and degasifier 
assembly, the mass of water in the tank, and all other apparatus within 
the tank insulation. The piston represents the water surface, which 
in the actual tank performs' work upon the steam during an in-surge. 
Therraodynamic equilibrium is assumed at all times. 

Equation //I is the well-known work equation: 



w s -Jp dv' 

where w - work (BTU/lbj, positive when work is done on the system 



P - system pressure (psia) 
v r specific volume of steam (, cubic feet/lb) 



Considering units, 




Equation #2 is the First Law of Thermodynamics: 



A e r q + w 



where 



or 



e = q +• w + e o 

e - specific internal energy (BTU/lb) 
q s heat flow (BTU/lb), positive when into the system 
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Equation //3 is the equation of state for a non-ideal gas: 



Pv = c ] _ 2 T 

where Z is defined as Pv_ (dimensionless compressibility factor.) 

C]T 

T - absolute temperature (°k) 

3 

C 1 r 1545 = .5957 lbg ft (gas constant) 

18.01*144 in lb 7 °R 

Equation #4 relates P to e and v empirically from steam table 
values (see Appendix I): 

P = 8.51 e + 52,590 v 2 - 30,700 v - 3161 



Equation #5 relates Z to P and v empirically, in the saturated 
and superheated region, from steam table values(see Appendix 1): 

Z = 1.715 x 10"^ P - 9.47 v 2 + 6.04S v - .5693 
Note: Equation #5 is not required mathematically, but is desirable 

for computer circuit simplicity. 



Equation #6 relates heat flow to temperature differences and is 
developed as follows: 

First, consider the sink as a system, 

d ^sink = m sink c sink dT s ink 
where '^sink = quantity of heat energy (BTU) 

m sink = mass of sinic 

T s i n k = temperature of sink (°R) 

c sirjk = s:i - rdc heat capacity (BTU/lb°R) 

Rearranging, 



dT 



sink 



d ^sink 



■^inkCgink 



6 



Integrating, 
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sink 


1 nk 




u sink c sink 


but 


*stea:i 


" ^sink = Q 


and 


(t ) 

v o' steam 


= ^ T °^sink s r o 


Then, 







T 

sink “ 



- 1 

m sink c sink 
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T 



o 



bow, considering the original system (the steam): 

- -K (T - Tsink) 
dt 

K represents the effective thermal conductance of the system boundary, 

and will be considered a constant in the equivalent system for any given 

surge. The units of K are (3TU/sec°Ii) . K is ordinarily the product 

hA, when each is determinable. 

Substituting for T 

d^ = - K (T + 4 - T 0 ) 

^ ;a sink c 3ink 



Since Q 



m s q , where m s = mass of steam (lb) 
da = - K (T - T c ) - K q 

^ m s ra sink c sink 



The problem which is now encountered is that of evaluating the 
equation constants. Since accurate design data are not available, these 
constants must necessarily be evaluated in an arbitrary and approximate 
manner. 

For the subject tank, can be taken as 3t>00 lb. This is 

approximated as follows: 
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mass of metal ~ 



2000 lb. 



mass of water 1000 lb. 

Total 3000 lb. 

The heat capacity of water at the pressure and temperature range 
involved is approximately 2.0 BTU/lb°R . The heat capacity of stainless 
steel (tank wall material) is approximately .13 BTU/lb°R . A weighted 
average heat capacity of the sink is computed as follows, 

c .. = 2 x 1000 + .13 x 2000 = .75 3TU/lb°R 

S 3000 

The mass of steam in the system will depend on the initial tank 
level, and saturation conditions, but rarely varies greatly from an 
average value of 175 lb. Properly, the mass of steam should be computed 
and Equation #6 modified for every run simulated, but for simplicity 
this average value is used for all surges. In all simulated surges, 
ra sink’ c sink> and m s are taken as constants, while the quantity K is 
adjusted to cause computed transient properties to coincide with experi- 
mental data. 

Equation //6 with constants evaluated becomes, 

d£ - - 5.72 x 10~ 3 K (T - T ) - .445 x 10~ 3 K q 

dt 

In Heaviside operator notation, the equation takes the form, 

q = - 1 £ 5.72 x 10" 3 K (T - T q ) + .445 x 10 _3 K qj 
? 

During all runs the assumption is made that the steam mass remains 
constant throughout the surge. Since the energy introduced as work 
during a surge is on the order of 3000 3TU, the maximum mass change 
would be in the neighborhood of 6 lb., if all of the energy were involved 
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in a phase change. Therefore, the percentage mass change can never be 
a significant value, especially since heat flows abundantly to the sink. 
Any phase change will significantly affect K however, because of the 
heat transfer mechanism involved (heat of vaporization). 

The six theoretical equations as developed and in their most useful 
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as follows : 
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3.07 x 10 4 


c 9 = .5693 


5.259 x 10 4 
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, -3 
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CHAPTrli 111 



APPLICATION OF i^UATIU.S TO THE ANALOG CU^UTta 



After the derivation of the six theoretical equations, a computer 
circuit must be designed to accomplish their simultaneous solution. In 
the circuit design .nany things must be considered; such as, scaling 
factors, follow-up time of multiplier components, amplifier drift, volt- 
age magnitudes, etc. A certain amount of latitude is allowed in circuit 
component selection, but definite limits are present due to stability 
considerations . 

The six basic equations in computer notation appear as below: 



n 

in 






itU 

//5 



w = - 1 [l.lll ( .IF dv dt ) | 

p >- dt 

e = .0333 q + -0333 w + e 

o 

T = 3.3574(~j 

P = 4-255 e ~ 5.1167 v + 4.3825 l.02v J - 52.6833 

7 = .5145 e + 3.024 V - 2.3675 ( .02 v 2 ) - 28.465 

“ = - i[. 1716 K T - .1716 K T + .000445 K q ] 

P L 



The bar over a quantity represents the quantity as a voltage. The 
detailed conversion of the theoretical equations to the above "machine " 
equations co. iprises Appendix II. 

Fig. 3 shows the resulting machine circuit in block diagram form. 
For the detailed circuit, see Appendix 11. Briefly, the circuit is 
designed to function as follows: 
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BLOCK DIAGRAM OF COMPUTER CIRCUIT 

FIGURE 3 




In the driving function circuit a cosine function is udaed elect- 
rically to a constant voltage, the resultant voltage transient very 
closely approximates the specific volume throughout a surge. The 
amplitude of the cosine function represents one-half the change in v , 
while the period of the function is twice the time of the surge (a t) . 

Within this same circuit the time derivative of specific volume, dv 

dt 

is also generated. 

A function multiplier is used to obtain the quantity , while a 

servo-multiplier produces the products P dv and P v. An integrating 

dt 

circuit integrates P dv as a function of time to produce specific work, 

— dt 

w . Another integrating circuit produces the time integrated specific 
heat flow, q , using temperature inputs . 

A summing amplifier sums e Q , w , and q to produce e. Another sums 
e, v, and v^ to produce P. A third sums P, v, and v^ to produce Z. 

A dividing circuit performs the division P v to produce T. 

Z 

The circuit required to solve the set of equations consists of 
eleven amplifiers, four sign-changing amplifiers, two function multipliers, 
one duo-channel servo-multiplier, and twenty-one 50,000 ohm potentio- 
meters. The net result is a rather complex circuit, but one which does 
not seem to be subject to further simplification without sacrificing 
necessary accuracy. Since instability in this type of computer is a direct 
function of the number of components used, circuit simplicity has been 
a constant objective. 

This phase, the designing and setting up of the circuit, requires 
a great deal of computer experience and proficiency. For the uninitiated, 
many trial and error situations are encountered, and with a circuit of 
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this complexity the process can be very time consuming. 

The following chapter wall deal with the problem of producing 
quantitive results , assuming the circuit has been designed, balanced 
and tested. 
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CHAPTER IV 



ANALOG COMPUTER SOLUTION OF THEORETICAL EQUATIONS 



After the analog computer circuit is designed, assembled, and rough 
qualitative results obtained, the next phase consists of producing 
accurate quantitative solutions. The first step in accomplishing this, 
consists of static adjustment. With various initial values of e Q , and 
v set, the e, P, Z and T computing circuits are adjusted to give steam 
table values for these quantities over the entire anticipated range. 

This step ofthe co puter set-up is very exacting and time consuming 
but is a necessary step prior to adjusting for dynamic accuracy. 

Prior to the dynamic adjustment, a hand solution of the system of 
theoretical equations is obtained for use as an adjustment reference. 

A set of initial conditions, and a driving function are arbitrarily 
selected, and using the Heun Method'' of numerical integration, a point- 
by-point solution of the equations is calculated (see Appendix III). 

The various thermodynamic quantities as calculated are plotted, the 
resulting curves serving as the final reference for dynamic adjustment 
(see Fig. 4) • 

Next, the computer is made to duplicate the hand solution. This 
is done by adjusting the w and q integrating circuits until, the w and q 
curves, as viewed on a properly calibrated Sanborn Recorder, are of the 
proper shape and magnitude (see Fig, 5). These adjustments in general 
are minor. 

The value of careful static adjustment is now apparent since no 
further adjustment of the e, P, Z and T circuits is found necessary. 



1 Special case of Runge-Kutta Method; see Kopal, Numerical Analysis, 
pp. 202-205, Wiley (1955; . 




>5 
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SANBORN RECORD gJG OF TEST RUN 



-p * 
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The curves of these quantities are non found to match the cori - esponain : 
hand solution curves very well (see Figs. 4 and 5 



Once the hand solution is duplicated, a reference for measuring K 
is now available. In the hand solution, the arbitrarily selected value 



readily observed (usually in the neighborhood of .100, varying slightly 



It is now possible to vary K over a large range of values by merely 
adjusting a single potentiometer, thus varying the integrated L curve, 
and all other curves of thermodynamic properties. 3y simulating the 
driving function for a given surge, K can be varied until the proper 
pressure curve is obtained (i.e. a pressure curve which matches the 
pressure curve resulting from an actual tank transient J. The next 
chapter discusses the investigation of actual surges by simulator, to 
determine values of K for these surges. Correlation of K with surge 
types is also discussed. 



of K is .9254. The corresponding setting of potentiometer a, is now 
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from day to day,). Since K and the a^ setting are proportional, K can 
be determined at any time by the simple relation, 



h = K 0 x (ayj 
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CHAPTLm. V 

IiXPERIldliNTAL itDSULTb AND CCNCLUolUJS 

Appendix IV is a compilation of data obtained from tne actual surge 
tank, eighteen positive surges are represented, spray being utilized 
in only one. It is evident that possession of many more runs would be 
desirable, but the present number must sui'fice since more could not be 
obtained. Fortunately these eighteen surges represent a lair cross- 
section of positive surges encountered. 

The procedure for simulating a given surge is outlined in Chapter 
IV and briefly is as follows: 

A proper driving function (specific volume} is set up on the comput- 
er, providing for the generation of a cosine curve of proper amplitude 
and period. Initial values of v, e, T, and a trial value of K are set, 
and a surge is generated. K is adjusted for several trials until the 
generated property curves match the corresponding actual tank curves 
(see Fig. 6 for Sanborn recording of this procedure). 

iiach of the eighteen surges is simulated in this manner, and the 
proper value of K -determined (see Table 1). kany attempts at correlat- 
ing K with various quantities (a v, A_v , iv , etc.) were fruitless, but 

v o At 

a fair correlation with surge duration <At} is possible (see Fig. 7). 

It is likely, in fact, that K is a function of several or all of the 
variables considered, but the effect of each seems small in comparison 
to the effect of At. 3y simple curve fitting, a correlating equation 
is derived, 

K = 3.75 + 2£_ 

t - 17 
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K iSTlLmNATICA FOR TYPICAL RUN 
20 
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Nothing is claimed for this equation except that it describes a K which 
aids in predicting pressure rises on the safe side (i.e. in most cases 
the predicted pressure rise wiil be somewhat larger than the actual 
pressure rise; . Column A of Table 1 lists the K for each surge as 
predicted by formula. 

In general, it can be observed that K does not vary greatly for 
any surge of over 30 seconds duration, and that it becomes nearly cons- 
tant for surges longer than 90 seconds. It is also true that pressure 
rise is not sensitive to moderate changes in K. The last colu.on of 
Table 1 lists errors encountered when a constant value of 3 - 75 , the 
limiting value of K, is used in the computer. 

The mechanis n of heat transfer for the actual tank can only be 
surmised. However, it is logical to expect greater turbulence for short- 
er, faster surges, resulting in somewhat larger values for K . This 
could be due to a larger heat transfer surface between the steam and 
water in the tank. The effectiveness of spray in increasing K is readily 
observed in run 1506, Table 1. V/ he re as a K of approximately 5-0 could 
have been expected without spray, the actual K observed is 13.8 . 

From the general characteristics of the K vs. Surge Duration curve 
(Fig. 7), it appears that K approaches a limiting value of about 3.75 
for surges of long duration. For geometrically similar tanks utilizing 
the same pressure range, the same limiting value of K could reasonably 
be expected. 

Realizing that a correlation exists between K and the driving 
function, it is now possible to predict pressure rises for all conceiv- 
able driving functions. This is done in the following mariner: 
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■i ijrror in predicting pressure rise (in psia^ v/hen using a constant 
value of K,( the Uniting value of 3-75;, rather than the exper- 
imental or formula values. 
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A driving function is e; tablis.ied, the proper value of K is set, 
a..d the resulting pressure rise recorded. In an attempt to use tiie 
principle of geometric similarity, dimensionless parameters are used 
when possible. Accordingly , a fa 'ily of cur/es is obtained relating 

a v , a p , and at (see Fig. 8; . 

v P. 

o . o 

Since the acoual tan!; surges vary fro. 27.5 seconds to 120 seconds, 
the predictions are only valid v/ithin tliis range . However it appears 
that safe extrapolations may be made . An attempt to justify limiting 
values of the prediction curves is based on the following reasoning: 

initially the system consists of ary saturated steam, «s a surge 
progresses 'work is done on the system, and heat will flow out provided 
a temperature difference exists. Examination of Fig. 9 shows that the 
isentx'opic compression process prescribes the practical upper limit of 
over-pressure. It is believed that a process following the saturation 
line prescribes the iower limit of over-pressure for this tank. Since 
the surge introduces colder water into the ta;ik, a potential sink exists 
which could conceivably cause the process to enter the wet steam region. 
However, a careful study of actual surges indicates that the stead 
remains superheated throughout an in-surge . Very long surges appear to 
parallel closely the saturated line with slight superheat. In addition, 
all surges appear to oe initially isentropic in nature . 

Based on the aoove reasoning and actual tesu results, it is believed 
that the prediction curves approach finite values as 1 1 approaches in- 
finity. Use limiting value of (for any given t v ) is that which 

r o V p 

would be obtained in a process following the saturation line. 
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TEMPERATURE - ENTROPY CHART (APPROX IM ATE) 

FIGURE 9 
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For very short surges, the li .itin~ cond^ti n appears to be o.io 

of isentro^ic co..ipression, thus establishing a -uximum value of,^. 

_ . : o 

for any given v . On the actual installation, curves of less than 10 
v D 

seconds duration are uid-i. ely quo to primary loop circulation ti.^e 
considerations. Thus the extrapolated curves for values of a t from 10 
seconds to 30 seconds would seem reasonable. Using the prediction 
curves (Fig. 8), Table 2 shows a comparison of predicted and actual 
pressure rises. 

It is believed that basic objectives of this project have been 
realised; however ;nany avenues for further study appear tc be open. 
Although a procedure for rational analysis has been indicated, it appears 
that a study of the effect of spray, and of negative surges, by these 
methods would be of value. 

In the field of heat transfer, a great deal remains to be done. 

The combining of all heat transfer mechanisms into a single coefficient 
is perhaps an over-simplification, although effective in this analysis. 

It is known that several thermal paths exist, and a more complete 
analysis could perhaps isolate the effects of each. Heat can flow from 
steam to metal, from steam to water (directly or by condensation;, and 
from metal to 'water. A comprehensive simulator study could provide 
knowledge in this relatively unexplored field. 
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rue eg. 


1 

; . .uJDicT^i, r 

(psiaj 


AC rUAL P 
(psiaj 


PREDICT 

-(2 


2302 


179 


174 


5 


2045 


317 


313 


4 


2128 


310 


3 j 0 


10 


1556 


236 


229 


7 


2229 


293 


303 


- 5 


1506 


302 


239 


13 


2023 


214 


242 


-28 


2010 


212 


232 


-20 


2243 


257 


261 


- 4 


2240 


130 


110 


20 


2247 


380 


314 


66 


2113 


150 


187 


-37 


1101 


139 


112 


27 


2255 


179 


152 


27 


1911 


420 


417 


3 


1658 


335 


310 


25 


I630 


329 


324 


5 


2039 

! 


362 


341 


21 
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APPENDIX 1 



DERIVATION OF EMPIRICAL RLLATluho FOR PuEJoURE 
aND C0EPIuIooI3ILlTY FACTOR 

oince e can be computed from the First Law, e = q + w + e Q 
and v is available as a driving function, it is convenient to express 
P as a function of e and v . This can bo done by constructing an 
empirical relationship from steam table values, which will De accurate 
within cex'tain limits. The pressure range which is considered signifi- 
cant for this problem is 1700 - 2700 psia, and no accuracy is claimed 
beyond these limits. 

Since e cannot be taken directly from the steam tables in the 
superheat region, a plot of enthalpy vs. pressure is made initially, 
showing lines of constant temperature and constant superheat (see Fig. 
10 ). 

Next, at various temperatures, various values of h , P , and v 

are listed. For each of these values e is coraputed from the relation 

e = h - 144 P v (see Table 3 )- From this taole, a new chart is 

J 

constructed relating e , P , v and T (see Fig. 11). Using this chart 
as a basis, initially an expression for e = f (P,v) is obtained in 
the following manner: 

The average slope of constant v lines is .1176. For any one line, 
e = ,1176 P +■ c-p.*. • The constant c-y is determined as a function of 

v (see Table 4 and Fig. 12;, c-jj = 1562 v + 538. The resulting 
expression for internal energy, e = .1176 P +- 1562 v +- 538 , when 
tested at the extreme limits is found to give considerable error. 
Although undesirable from the standpoint of computer circuitry, an 
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FI6URE 10 




TABLji 3 



T (° F ) 


h ( 3 TU/lb) 


P (psia) 


v ( ft^/lbj 


P V 


e ( 3 rU / lb ) 


630 


1136.1 


1700 


.2553 


8 o .5 


1105.6 




1167.7 


i 860 


.2296 


76.6 


1091.1 




1145.8 


1900 


.2040 


71.6 


1674-0 


(sat . ) 


1141.1 


1919.3 


.1992 


70.8 


1070.3 


b! r J 


1185.7 


1800 


.2407 


80.5 


1105.4 




1167.0 


I 960 


.2168 


76.4 


10 / 0.6 




1145.6 


2000 


. 1/36 


71.7 


1073.9 


(sat.) 


1130.6 


2059.7 


.1798 


68.7 


1061.8 


650 


1167 . y 


2000 


.2058 


76.2 


1090.8 




1146.3 


2100 


.1846 


71.8 


1074.5 




1121 . U 


2200 


.1633 


06.4 


1054.6 


(sat . ) 


1118.5 


2208.2 


.1616 


66.2 


1052.3 


660 


1167.7 


'2100 


.1962 


7 u .3 


1091.4 




1147.3 


2200 


.1763 


72.2 


1075.6 




1123.8 


2300 


.1575 


67.3 


1056.5 


(sat. ) 


1104.4 


2365-4 


.1442 


63.O 


1041.4 


670 


1150.0 


2300 


.1702 


72.5 


1077.5 




1123.2 


2400 


.1526 


67.8 


1060.4 




1099.8 


2500 


.1342 


62.2 


1037.6 


(sat . ) 


1087.7 


2531.8 


.1277 


59.7 


1028.0 


630 


1132.3 


2500 


.1484 


68.7 


1063.6 




1107.1 


2600 


.1319 


63.5 


1043.6 




1072.3 


2700 


.1137 


56.8 


1016.0 


(sat. ) 


1067.2 


2703.1 


.1115 


56.0 


1011.2 


690 


1156.6 


2500 


.1594 


73 . S 


1082.8 




H 37.3 


2600 


• lt 47 


69 • 8 


1067-5 




1114.3 


2700 


.1299 


6/ r .8 


lu 49. 3 


700 


1160.6 


2600 


.1549 


74.6 


1086.0 




1142.5 


2700 


.1415 


70.7 


1071.8 
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TaJlL 4 



V 


e 


P 


.1176 p 


C 14 


.12 


1021.0 


2620 


308.5 


712.5 


.13 


1030.4 


2510 


295.5 


735.0 


.14 


1033.2 


2409 


284.0 


754.2 


.15 


1045.4 


2312 


272.0 


773.4 


.16 


1051.5 


2217 


261.0 


790.5 


.17 


1056.8 


2140 


252.0 


804.3 


.IS 


1062.0 


2057 


242.0 


820.0 


.19 


1066.5 


1982 


233.5 


333.0 


.20 


1070.5 


1915 


225.5 


845.0 
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FIGURE 12 






additional modification is necessary, making e = g( P, v , v 2 ). 

n 

To accomplish this, the expression P - 6. 51 e - c^ v + c. 
+ c-^y is set up. Substituting actual steam table values at three 
points ( P = 1900, 2100, 2300) three simultaneous equations in c-^ , 
c-^. , and c-jy, result. Solving for these constants and substituting, 
the final equation for P is : 

P - 8.51 e + 52,590 v 2 - 30,700 v - 3161 



The calculation of an empirical relationship, L - f (p, v) follows 

similar lines. First Z vs. P is plotted from the steam tables, 

showing lines of constant v (see Fig. 13). The average 3lope of the 

4 

constant v lines is 1.715 x 10“ , the equation for any one line being 
Z - 1.715 x 10~ 4 P + c-^.j, . c-^ is determined as a function of 

v , c . = 2.943 v - .31S (see Table 5 and Fig. 14) . The expression 
for Z as a function of P and v is, Z - 1.715 x 10 P +■ 2.943 v 
- .318 

p 

Once again a v term is found necessary for desired accuracy, and 
using the same procedure as before the resulting expression for Z is: 

Z = 1.715 x 10 -4 P - 9.47 v 2 + 6.048 v - .5693 
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TA3LE 5 



I 

V 


Z 


P 


-4 

1.715 x 10 P 


C 18 


.12 


. 464O 


2613 


.4475 


.0165 


.13 


.4845 


2509 


.4300 


.0545 


.14 


.5030 


2410 


.4130 


.0900 


.15 


.5225 


2312 


.3960 


.1265 


.16 


.5380 


2211 


.3790 


.1590 


.17 


.5510 


2140 


.3670 


.1840 


.18 


.5640 


2060 


.3530 


.2110 


.19 


.5770 


1985 


.3405 


.2365 


.20 


.5390 


1915 


.3280 


.2610 


.21 


.6005 


1850 


.3170 


.2835 



t 



3S 
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APPENDIX II 



DliiillVA CF Mj-tChlKK il ^U n To. OKS 



Scaling factors are assigned as follows: 



£rr.ax 

* max 



3000 

50 



60 



o< 



_ V 



max - 



max 



50 



.01 



= T nax = 1500 

= 50 

x max 



30 



-x _ = e 



nax 



max 



1500 

50 



30 



* = z i 



'max 



Z. 



max 



50 



.02 



<x„ = w, 



w 



max 



w, 



max 



^O 

50 



o< , 



maax 

%iax 



^2 

50 



o( 



t = - 



1 

1 



Conversion of Theoretical Equations to Machine Equations 

» . - I [e 12 P (g) * J 



w = 



VI - 



- 1 \ .1851 (.1 P dv| dt] 10 

p L dt 1 j — 



-if 1.111 (.1 P dv\ 

p L V dt / 



at ] 



40 



ff2. 



e - q +• w + e c 



e = ^ q q a ^t L W a. <^e e 0 






Of , 



c*, 



= .0333 q + .0333 w + e c 



//3 



T m P v 
CqZ 



#4 



#5 



/.01]T 


■ \ .01 P V 


\ .02/ 


ci«* z <*r \__.02 z 


7 - 


3.3574 i .01 P_V ] 




1 .02 L J 


P = 


C2 e - v + c 4 A 


P = 


c<e c 0 e “ °^v c _2_ v 






P = 


4.255 e - 5.1167 ~ 


Z = 


C£j P +■ C7 V - Cg V 


z = 


C^O<p P -p- c^c< v V 




c< z exTz, 


z = 


,5145 P + 3.024 V 



- c. 



+ 50 2jgvffv ( * u2 v ) ■ S 

<*p ^ 



c 9 



- 50 c R oP v o< v (.02 V J - £ 9 . 

7 . 



C<z 



- 2.3675 (.02 V 2 ) - 28.465 



// 6 



q = - 



q = - 



lfc i0 K (T - T 0 ) + C U K q_ 

P ^ 

C 10 K ^ ^ ' V -4- c 1 1. K 



0< 






C< t 
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q = - 1 .1716 K r - .1716 K r 0 + .O004A5 K q 

P U 



Calculated potentiometer settings, and components selected: 



a l 


= 


D. F. 


dependent 


a 8 


= 


.4255 


a 15 


= 


.4383 


a 2 


= 


ii 


)t 


a 9 


= 


.6864 


a l6 


- 


.6864 


a 3 


a 


it 


tl 


a iu 


= 


. 5H7 


a 17 


= 


.5000 


a 4 


= 


.1111 




a ll 




.6333 


a l8 


= 


.Ou89 


a 5 


- 


.3333 




a 12 


s 


.3024 


a 19 


r 


Iq. 

45 


a 6 


— 


.3357 




a 13 


= 


.5145 


S 20 


= 


.2368 


a? 


= 


.3333 




a 14 


= 


K 


a 21 


- 


.1182 



R = 


1 


M 


0 

Hj 

1-* 

II 


. lytf 


R 45 “ 


10 


11 


O 

I- 1 * 

& 

R 


2 j£ 


a 46 = 


10 


M 


C fS = 


• l/'f 


R f = 


1 


Li 


C flO* 


.l/*f 


II 

cm 


.1 


M 


c f28“ 


c 

c 

'S 



Rf 30 = lu M 

a 60 * *1 K 
R 5 = .1 M 
R^ - 10 M 

a 68 “ * 1 
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S, , S, ,s 5 ,$4 SIGN CHANGERS 

COMPUTER SCHEMATIC DIAGRAM 

FIGURE. I JT 
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APPENDIX 111 



HAND SOLUTIuh OF THEORETICAL EQUATIONS 



The first step in obtaining a hand solution consists of combining 
the six theoretical equations to form two dependent differential 
equations : 

Substituting from equation #2 into equation //4, and then into 
equation // 1 one obtains, 





<dt>7 

d-t- " 


C-L 


Since 


V = 


v'«,. f «« “>* 


then 


dv . 

aa - 


” U 0 to Slid tut 


and 


du7_ 





- ( -V 1 dy 
c t J dt: 



d t 



-f 1 17 



_ C 



r?V + 






Cj. J 



Substituting from equations // 4 and //5 into equation #3, the 
resulting expression for T becomes, 

V Ccj + v/N e 0 - j 



T- 



C, 



^6 



c* 






‘C.C* 



ia 

c, c*c 6 



It is now convenient to assign values to the driving function, 



initial conditions, and to evaluate the constants: 



u o = 


.019 ft 3 /lb 


Oj = 


.0418 rad/sec 


V 

mean • 


= .172 ft 3 /lb 


e o “ 


1066 3TU/lb 


Cl = 


.5957 



c 2 = 8.51 

c-} - 3.o7 x lu 4 

c 2 = 5.259 x 10 4 

c^ = 31 ol 

c 6 = 1.714 X 1U~ 4 



Cy ** 6.U48 
Cq =9.47 

Cy = .5593 

c 12 = * 18 53 
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Evaluating, 

dt-j 
d t* 



.2b3 *io sm .o4io^ [<j -i-"u7 - 3faiO v + ioioO v ■+ 6 94J 1 



-j- _ ?8oo v [] <^ -ruj - 3 6 to v -t- £>i0o v 1 t 694] 
[(jtw t 540 V - 3/0 V~ ■+ 304] 

Equation ij 6 is, 

dt = ~ C '“ KT + C » k:T o - <4 K % 



Arbitrarily assigning a value to h , evaluating constants, and 
substituting the already deriveu expression lor T into equation // 6, 




K = 




.9254 

5.72 x 1CT 3 
4.45 x 10~ 4 
1094 • 5 °a 



-51.8 7 “36/0 V t-btSoV 2 * 694] 

r Xtf + 54 O V * 310 v L 1 304 1 



4-/2*/o'(g t 5.7 7 



The two resulting dependent i'irst order differential equations 
become , 



r AT, r^+W+fcrJ -4,/Zno^ 1 5.110 
dl C<3 ♦ W + < 3 ] 

Lr w + ^J 



where K, } K z , A], are tabulated values of V , the driving function 
(see Table 6). 
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TABLE 6 



Time 


U 0 COS<ut 


V 


K 1 


K 2 


k 3 


*4 


0 


.0190 


.1910 


-9.90 


230.5 


395.7 


0 




5 


,0l86 


.1906 


-9.88 


2^1.0 


395.5 


2.6l 


x 10~^ 


10 


.0174 


.1893 


-9.82 


232.0 


395.1 


5.10 


11 


15 


.0154 


.1873 


-9.71 


2^5.0 


394-3 


7.38 


11 


20 


,ui27 


.x847 


-9.58 


239. U 


393.1 


9.33 


ft 


25 


.0095 


.1315 


-9.41 


242.5 


391.8 


10.87 


M 


30 


.0059 


.1779 


-9.24 


247.5 


390.3 


11.90 


tl 


35 


.U02u 


.i74o 


-9.02 


253.0 


388.6 


12.47 


1! 


• 40 


-.0020 


.1700 


-6.82 


259.5 


336.6 


12.47 


It 


45 


-.0059 


.1661 


-8.62 


264.5 


335.2 


il.9o 


!t 


5o 


- . u095 


.1615 


-3 . 57 


272.3 


383.1 


10.87 


n 


55 


-.0127 


.1593 


-8.26 


276.0 


382.1 


9.53 


11 


60 


-.0154 


.1567 


-8. 13 


261.0 


381.0 


7.38 


11 


65 


- . oi74 


.1547 


-0 . U2 


283.7 


380.1 


5.10 


11 


70 


-.0186 


.1534 


-7.95 


285.0 


379.6 


2.61 


n 


75 


-.0190 


.1550 


-7.94 


236.5 


379.5 


0 




so 


-.0166' 


.1534 


-7.95 


289 .0 


379.6 


-2.61 


11 


65 


-.0174 


.1547 


-8.02 


283.7 


380.1 


-5.10 


11 


90 


-.0154 


.1567 

■ 


-3.13 


281.0 


381.0 


-7-38 


it 


95 


-.0127 


.1593 


-8.26 


276.0 


582.1 


-9.33 


11 


100 


-.0095 


.1615 


-8. 37 


272.3 


333.1 


-10.87 


11 
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The next step consists of numerically integrating the above 
equations to obtain point-by- oint values of ^ ^ , \tj , and dp. 

The Heun L'ethod of numerical integration is used, the general form 



being: 



where 



where 



% ' 


f (vj 




,v) 




TV “ 


3 (w 


> l Z 


y) 




<*> 

% 


V' 


+ 


\ ( 9,..-, 


, (h-ii 
1 ‘ ) 


m 










~~ 




i 


11 4- 










t / . 




Wc 


w c-, 


■+ - 


5 ^ “k. 


+ ) 


(A 

07. = 




+ 


h ( W. \ 



h is the interval between points and is selected as 5 seconds. 
Table 7 lists the resulting values for Q. , , \<J , and xij 

obtained over a time of 1U0 seconds (21 points). 



Next, the resulting values of and Tj are substituted back into 
the original theoretical equations, resulting in point -by-point values 
for the various properties throughout the run. These values are tab- 
ulated (see Table 8), ana are plotted as a function of tine in Fig. 4, 
Chapter IV. 
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TABLE 7 



HAND SOLUTION (K =.926 5/ 



i 


t 


q i 




W 

Jl 


w. 

1 


0 


0 


0 


0 


0 


0 


1 


5 


-.0273 


-.Olio 


.1509 


.0604 


2 


10 


- . 1083 


-.0215 


.6011 


.1197 


3 


15 


-.3410 


-.0716 


1.3431 


.1772 


4 


20 


-.3431 


-.1292 


2.3614 


.2302 


5 


25 


-1.5479 


-.1527 


3.6260 


.2757 


6 


30 


-2.4188 


-.1957 


5.0948 


.3119 


7 


35 


-3.4896 


-.2327 


6.7173 


.3372 


8 


40 


-4.7612 


-.2760 


8.4311 


.3483 


9 


45 


-6.1402 


-.2756 


10.1536 


.3407 


10 


50 


-7.5859 


-.3026 


11.8070 


.3206 


11 


55 


-9.1073 


-.3062 


13.3039 


.2781 


12 


60 


-10.6548 


-.3126 


14.5556 


.2226 


13 


65 


-12.1476 


-.2845 


15.4969 


.1539 


14 


70 


-13.4633 


-.2438 


16.J765 


.0779 


15 


75 


-14.6571 


-.2317 


16.2714 


0 


16 


SO 


-15.6849 


-.1794 


i6.0843 


-.0748 


17 


85 


-16.5174 


-.1536 


15.5423 


-.1420 


18 


90 


-17.1731 


-.1087 


14.6944 


-.1972 


19 


95 


-17.5063 


-.0246 


13. 6080 


-.2374 


20 


100 


-17.4614 


.0426 


12.3475 


-.2659 



AS 



TABLE 8 



PROPERTY VALUES FRULI HAND SOLUTION 



t 


e 


V 


p 


2 


T 


0 


1066. U 


.1910 


- 

1966.9 


— — — | 

.5774 


1092.1 


5 


1066.12 


.1906 


1969.7 


.5774 


1091.3 


10 


1066.49 


.1893 


1986.4 


.5773 


1093.4 


15 


1067.0 


.1873 


2015.4 


.5766 


1098.9 


20 


1067.52 


.1847 


2047.0 


.5760 


1101.9 


25 


1068.08 


.1815 


2086.9 


.5747 


1106.5 


30 


1063.68 


.1779 


2134.2 


.5734 


1111.4 


35 


106V.23 


.1740 


2190.1 


.5718 


1118.8 


40 


1069.67 


.1700 


2243.0 


.5699 


1123.1 


45 


1070.01 


.1661 


2297.3 


.5679 


1127.9 


50 


1070.22 


.1615 


2361 . 5 


.5653 


H32. 6 


55 


1070.20 


.1593 


2391.9 


.5638 


1134.3 


60 


1069.90 


.1567 


2427.1 


.5616 


1137.0 


65 


1069.35 


.1547 


2447.0 


.5597 


1135.4 


70 


1068.61 


.1534 


2459.6 


.5578 


1135.4 


75 


1067.61 


.1530 


2458.1 


.5560 


1135.5 


SO 


1066.40 


.1534 


2440.8 


.5546 


1133.2 


85 


1065.02 


.1547 


2390.1 


.5499 


1128.6 


90 


1063.29 


.1567 


2370.9 


.5520 


1129.9 


95 


1062.10 


.1593 


2323.0 


.5520 


1125.5 


100 


1060.89 


.1615 


2282.1 


.5517 


1121.6 
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APPENDIX IV 



ACTUAL TANK TRANSIENT DATA 



Run #2247 
Time 


mo 

Pressure 


=, 169.06 lb. 

Level 


Volume 


Spec. Vol. 


0 


2002 


32.2 


31.7 


.1375 


5 


2012 


32.4 


31.6 


.1867 


10 


2042 


33.2 


31.2 


.1343 


15 


2094 


34.6 


30.4 


.1797 


20 


2168 


36.7 


29.1 


.1720 


25 


2233 


39.7 


27.3 


.1612 


30 


2280 


41.8 


26.1 


.1542 


35 


2305 


43.1 


25.4 


.1900 


40 


2316 


43.7 


25.1 


.1485 


Run #2240 


m 

0 


= 167.46 lb. 






0 


2002 


32.8 


31.4 


.1875 


5 


2007 


33.0 


31.3 


.1870 


10 


2026 


33.5 


31.0 


.1850 


15 


2047 


34.1 


30.6 


.1827 


20 


2071 


34.9 


30.2 


.1803 


25 


2091 


35.6 


29.7 


.1773 


30 


2104 


36.3 


2y . 5 


.1760 


35 


2112 


36.8 


29.0 


.1731 


40 


2112 


37.0 


28.9 


.1726 



50 



Run # 2229 



167.13 lb. 



m 



o 



Time 


Pressure 


Level 


Volume 


Spec . Vol . 


0 


1992 


32.5 


31.5 


.1885 


10 


2007 


33.0 


31.3 


.1873 


20 


2033 


33.9 


30.8 


.1343 


30 


2070 


35.2 


30.1 


.1800 


40 


2117 


36.8 


29.0 


.1735 


50 


2170 


38.6 


23.0 


.1675 


60 


2220 


40.4 


27.0 


.1615 


70 


2264 


41.8 


26.1 


.1560 


80 


2295 


42.7 


25.7 


.1537 



Run # 2128 m o = 169.05 lb. 



0 


2000 


1 

32.3 


31.7 


.1875 


10 


2012 


32.7 


31.5 


.1861 


20 


2037 


33.4 


31.0 


.1832 


30 


2072 


34.5 


30.4 


.1796 


40 


2114 


36.0 


29.5 


.1743 


50 


2162 


37.8 


28.5 


.1684 


60 


2211 


39.5 


27.4 


.1620 


70 


224 4 


41.0 


26.6 


.1572 


30 


2288 


42.0 


26.0 


.1537 


90 


2299 


42.6 


25.7 


.1520 


100 


2299 


42.7 


25.6 


.1514 

1 
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Run #ai3 
Time 


m o = 
Pressure 


172.7 lb. 
Level 


Volume 


bpec.Vol. 


0 


2000 


31.2 


32.4 


.1878 


2.5 


2002 


31.4 


32.3 


.1870 


5 


2019 


31.8 


32.0 


.1854 


7.5 


2054 


32.4 


31.6 


.1830 


10 


2090 


33.3 


31.0 


.1795 


12.5 


2124 


34.3 


3o.5 


.1765 


15 


2155 


35.3 


29.9 


.1730 


17.5 


2172 


36.2 


29.4 


.1703 


20 


21 82 


36.9 


29.0 


.1630 


22.5 


2136 


37.3 


28.7 


.1662 


25 


2187 


37.6 


28.5 


.1650 


27.5 


2187 


37.9 


23.4 


.1645 


30 


2187 


33.0 


28.3 


.1640 


Run #1101 


m = 
o 


169.5 lb. 






0 


1990 


31.6 


32.1 


.1892 


2.5 


1995 


31.8 


32.0 


.1888 


5 


2009 


32.3 


31.7 


.1870 


7.5 


2039 


32.8 


31.4 


.1854 


10 


2066 


33.5 


31.0 


.1830 


12.5 


2085 


34.4 


30.4 


.1793 


15 


2097 


34.9 


30.1 


.1775 


17.5 


2102 


35.3 


29.9 


.1765 


20 


2102 


35.5 


29.8 


.1760 



52 



Run it 2255 



175.0 lb. 



Time 


Pressure 


.Level 


Volume 


Jpec . Vol . 


0 


1995 


30. 0 


33.0 


.1385 


2.5 


2002 


30.1 


33.0 


.1885 


5 


2011 


3o.3 


32.9 


.1880 


7.5 


2032 


30.6 


32.7 


.1870 


10 


2043 


31.1 


32.4 


.1850 


12.5 


2055 


31.7 


32.0 


.1828 


15 


2094 


32.5 


31.6 


.1805 


17.5 


2121 


33.3 


31.1 


.1775 


20 


2134 


34.0 


30.7 


.1755 


22.5 


2143 


34.6 


30.3 


.1730 


25 


2146 


35.0 


30.0 


.1715 


27.5 


2147 


35.3 


29.9 


.1710 


30 


2147 


35.5 


29.8 


.1703 


Run #1506 


m o 


= 167.7 lb. 






0 


1952 


30.8 


32 . 6 


.1944 


5 


1958 


3^.9 


32.5 


.1940 


10 


1988 


31.7 


32.0 


.1908 


15 


2032 


33.0 


31.3 


.1867 


20 


2082 


34.4 


30.4 


.1812 


25 


2130 


35.7 


29.7 


.1770 


30 


2174 


37.1 


28.8 


.1713 


35 


2217 


33.5 


28.0 


.1670 


40 


2238 


39.7 


27.3 


.1629 


45 


2241 


40.0 


27.1 


.1616 
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Run #1556 mo = 159.7 lb. 



Time 


Pressure 


Level 


Volume 


Spec . Vol . 


0 


1393 


31.2 


32.4 


.2030 


5 


1394 


31.3 


32.3 


.2022 


10 


1923 


32.0 


31.9 


.2000 


15 


I960 


33.0 


31.3 


.1962 


20 


1998 


34.2 


30.5 


.1912 


25 


2034 


35.4 


29.3 


.1868 


30 


2063 


36.5 


29.3 


.1837 


35 


2086 


37.3 


28.7 


.1798 


40 


2101 


38.0 


28.3 


.1774 


45 


2114 


33.7 


28.1 


.1761 


50 


2122 


39.1 


27.5 


.1723 


55 


2122 


39.4 


27.4 


.1718 


Run //201C 
0 


ra Q 

1930 


- 164.6 lb. 

31.0 


32.5 


.1975 


5 


1933 


31.0 


32.5 


.1975 


10 


1942 


31.1 


32.5 


.1975 


15 


1962 


31.5 


32.1 


.1951 


20 


1938 


32.3 


31-7 


.1927 


25 


2021 


33.2 


31.2 


.1898 


30 


2057 


34.4 


30.5 


.1855 


35 


2094 ' 


35.4 


29.8 


.1012 


40 


2124 


36.4 


29.3 


.1782 


45 


2146 


37.0 


28.9 


.1757 


50 


2159 


37.6 


28.6 


.1740 


55 


2162 


37.9 


28.4 


.1727 



54 



Run #2023 



ru 0 — 166.6 lb. 



Tune 


Pressure 


Level 


Volume 


Jpec .Vol. 


0 


1946 


31.0 


32.5 


.1952 


5 


1949 


31.0 


32.5 


.1552 


10 


1962 


31.2 


32.3 


.1940 


15 


1982 


31.7 


32.1 


.1923 


20 


2008 


32.5 


31.6 


.1899 


25 


2040 


33.4 


31.0 


.1862 


30 


2074 


34.5 


30.4 


.1827 


35 


2111 


35.7 


29.7 


.1784 


40 


2144 


36.5 


29.2 


.1753 


45 


2168 


37.2 


28.8 


.1730 


50 


2180 


37.7 


28.5 


.1711 


55 


2186 


38.1 


28.3 


.1700 


60 


2188 


38.3 


28.2 


.1694 



Run #1630 m - 190.3 lb. 



0 


1938 


20.0 


39.0 


.1967 


10 


1952 


20.6 


38.6 


.1948 


20 


1982 


21.5 


38.1 


.1920 


30 


2022 


22.8 


37.3 


.1880 


40 


2074 


24.8 


36.1 


.1830 


50 


2132 


26.7 


35.0 


.1775 


60 


2184 


27.8 


34.4 


.1722 


70 


2224 


29.8 


33.2 


.1670 


GO 


2252 


31.4 


32.2 


.1623 


90 


2262 


32.6 


31.5 


.1589 



55 



:iun //1911 



134.7 lb 



Time 


Pressure 


Level 


Volume 


opec . Vol . 


0 


lyio 


23.2 


37.1 


.2007 


2.5 


1920 


23.6 


36.8 


.2004 


5 


1935 


23.8 


36.7 


.1993 


7.5 


1958 


24.2 


36.6 


.1974 


10 


1990 


25.4 


35.8 


.1950 


12.5 


2027 


26.0 


35.4 


.1914 


15 


2072 


27.5 


34.5 


.1872 


17.5 


2117 


28.8 


33.8 


.1830 


20 


2161 


30.0 


33.1 


.1788 


22.5 


2207 


31.8 


32.0 


.1737 


25 


2252 


33-0 


31.3 


.1690 


27.5 


2287 


34.5 


30.4 


.1645 


30 


2312 


35.7 


2V.7 


.1607 


32.5 


2323 


36.7 


29.1 


.1572 


35 


2327 


37.3 


28.7 


.1553 


37.5 


2327 


37.5 

' 


28.6 


.1548 



56 






Run , 1*1658 n 0 - 184.1 lb. 



Tine 


Pressure 


Level 


Volume 


Opec . Vul . 


0 


1972 


26.2 


35.3 


.1918 


10 


1972 


26.4 


35.2 


.1916 


20 


1982 


26.5 


35.1 


.1906 


30 


2002 


27.4 


34.5 


.1877 


40 


2045 


28.9 


33.7 


.1330 


50 


2092 


30.3 


32.6 


.1775 


60 


2147 


32.5 


31.6 


.1713 


70 


2200 


34.6 


30.4 


.1651 


80 


2249 


36.3 


29.3 


.1591 


90 


2278 


37.7 


23.5 


.1544 


95 


2282 


33.2 


2S.2 


.1532 
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Run //22 43 



158.3 lb. 



Tine 


Pressure 


Level 


Voluae 


Spec . Vol. 


0 


1903 


31.8 


31.9 


.2017 


2.5 


1903 


31.8 


31.9 


.2017 


5 


1907 


31. 8 


31.9 


.2017 


7.5 


1912 


32.0 


31.8 


.2008 


10 


1923 


32.3 


31.7 


.2000 


12.5 


1948 


32.3 


31.4 


.1982 


15 


1987 


33.6 


31.0 


.1958 


17.5 


2033 


34.6 


30.3 


.1912 


20 


2QS4 


35.8 


29.6 


.1869 


22.5 


2125 


37.0 


28.9 


.1825 


25 


2150 


38.0 


23.3 


.1787 


27.5 


2162 


38.6 


27.8 


.1755 


30 


2164 


39.2 


27.6 


.1742 


32.5 


2164 


39.4 


27.5 


.1735 


35 


2164 


39.5 


27.4 


.1730 
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Run #2302 ra Q - 162.5 lb. 



Time 


Pressure 


Level 


Volume 


Spec.Vol. 


0 


1883 


30.0 


33.0 


.2032 


10 


1891 


30.2 


32.9 


.2025 


20 


1902 


30.5 


32.8 


.2020 


30 


1918 


31.0 


32.5 


.2000 


40 


1935 


31.6 


32.1 


.1975 


50 


1955 


32.3 


31.7 


.1>52 


60 


1975 


33-1 


31.2 


.1921 


70 


1996 


34.0 


30.7 


.1890 


SO 


2017 


34.3 


30.2 


.I860 


90 


2033 


35.5 


29.8 


.1835 


100 


2054 


36.1 


29.4 


.1810 


110 


2062 


36.4 


29.3 


.1804 


120 


2062 


36.5 


29.2 


.1798 



59 






ilun // 2045 = 167.0 lb. 



Time 


Pressure 

— 


Level 


Volu,-ie 


Spec. Vol. 


0 


1902 


32.0 


31.8 


.190A 


5 


1990 


32.5 


31.6 


.189 A 


10 


2005 


32.8 


31. A 


.1880 


15 


2025 


33.5 


31.0 


.1856 


20 


205A 


3A.5 


30. A 


.1823 


25 


2091 


35.5 


29.8 


.1785 


30 


2130 


36.6 


29.1 


.17A6 


35 


2166 


37.5 


23.6 


.1710 


AO 


2197 


36.5 


28.0 


•167A 


A 5 


2227 


39.5 


27. A 


.16A0 


50 


2252 


AO. 5 


26.9 


.1610 


55 


2275 


a. 3 


26.5 


.158A 


6U 


2290 


A1.9 


26.1 


.1562 


65 


2295 


A2.1 


26.0 


• 15A8 


70 


2295 


A2.5 


25.7 


' .1539 



60 



Hun ;/2039 
I i :ne 


"0 

Pressure 


= 165.5 lb. 
Level 


7olu;:ie 


Jpec.7ol. 


0 


1981 


32.7 


3i.5 


.1905 


2.5 


1991 


32.' 


31.4 


.1902 


5 


1996 


33.0 


31.3 


.1896 


7.5 


1995 


33.2 


31.2 


.1888 


10 


2003 


33.6 


31.0 


.1874 


12.5 


2025 


36.1 


3o.6 


.1854 


15 


20/+5 


36-7 


3o.3 


.1632 


17.5 


2063 


35.4 


29.6 


.x802 


20 


2093 


36.1 


27.4 


.1770 


22.5 


2137 


37.8 


23 . 5 


.1734 


25 


2177 


33.2 


26.2 


.1698 


27.5 


2217 


3. '.2 


27.6 


.1662 


30 


2252 


40.3 


27.0 


.1626 


32.5 


2286 


41.2 


26.5 


.1590 


35 


2312 


42 . 4 


2>.8 


.1560 


37.5 


2322 


43 • 1 


2>.4 


.1540 
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